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A recent measurement by the DO collaboration finds a like-sign di-muon charge asymmetry in the 
B system that is roughly 3a larger than the value predicated by the Standard Model. This suggests 
new physics contributing to B — B mixing. For the current central value of the CP asymmetry, 
the required size of TI2 is larger than Standard Model estimates of this quantity. In this paper, 
we will explore the constraints on new physics contributions to ri2. We show that there are two 
dimension six operators of Standard Model fields in the electroweak Hamiltonian whose coefficients 
are not constrained enough to rule out possible contributions from new physics. We argue that a 
more precise measurement of T{Bs)/T{Bd), which is possible with currently available data, could 
either support or strongly constrain the existence of new physics in ri2. 



Both the DO and the CDF collaborations have mea- 
sured the like-sign di-muon charge asymmetry in the B 
system 
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Using 1.6 fb-i of data CDF obtained [T] A^^ = (8.0 ± 
9.0 ± 6.8) X lO"'^, while the DO collaboration recently 
reported [2] a resuh of A^^i = (-9.57 ± 2.51 ± 1.46) x lO"^ 
with 6.1 fb~l of data. Combining these results, one finds 
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.5 ±2.8) X 10" 
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This value is about 3(t away from the Standard Model 
(SM) prediction of A^^ = -0.2 x lO^^. 

Since these measurements are blind as to which flavor 
of B meson produced the two muons, receives contri- 
butions from the scmileptonic CP asymmetries of both 



_B, and B^ mesons, which we will call 



and 



re- 



spectively. The relation between and the afj is given 

by m 



Ali = (0.506 ± 0.043)afi + (0.494 ± 0.043)a^i , 
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DO [3] has also measured the semileptonic CP asymmetry 
a|] directly, albeit with large uncertainties 



(-1.7±9.1) X 10" 
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One can convert the di-muon charge asymmetry into a 
measurement of the semileptonic asymmetry of the Bg 
system using input from the i?^ system. If one assumes 
no new physics contribution to B^ mixing, one finds 
(combining with the explicit measurements) 



Ki)sM.^ =-(12-2±4.9)xl0 
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while using the measurement a^j = (— 4.7±4.6) x 10 
one finds 



Ki)„<^,_= = -(9-2±4.9) 

si 



X 10" 



While it is probably too early to tell if this discrepancy 
is due to physics beyond the Standard Model or fluctua- 
tions in the data, it is certainly interesting to understand 
what the implications of this measurement are for new 
physics (NP). 

There are two amplitudes each that characterize mix- 
ing in the Bq systems (g — s,d): the off-diagonal element 
of the mass matrix and the off-diagonal element of 
the decay matrix P^j- Only the relative phase 0'^ be- 
tween these two amplitudes is observable, such that one 
can choose the three real parameters |-/Wj'2|, \^i2\ a-nd (j)'^ 
to describe the physics. In terms of these parameters, 
the semileptonic asymmetry is given by 
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What are the values of the three parameters iMfjl, 
jPfal and (j)'' in the Standard Model? This question is 
not easy to answer, since both \Ml2 \ and jr'jl depend on 
non-perturbative physics, which is notoriously difficult 
to determine. However, much effort has been devoted 
to this problem, including lattice calculations to deter- 
mine required bag parameters. We use the calculations 
from [5|, supplemented with updated values for the de- 
cay constants and the bag parameters, obtained from [3]. 
Adding the uncertainties quoted in these to references in 
quadrature, we find 
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(9.8± l.l)ps-i 
(0.049 ±0.012)ps"^ 
(0.04 ±0.01) . 



(8) 



It should be noted that the calculation of Ff 2 relies on an 
operator product expansion, even though the energy re- 
leased is only mf, — 2mc ~ 2 GeV. Thus one might be wor- 
ried of the convergence of the expansion performed [7] . 

These parameters can also be related to other physical 
observables. In particular, they determine the mass and 
width difference between Bg and Bs mesons, as well as 
(6) the time- dependent CP asymmetry S^^. The relations 
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where we have assumed 



aTg[-Vt,V;jV,sV:,] « 
servables one finds 
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and 



0. In terms of these three ob- 
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The measured values for these three observables are 

[in]i2 
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= (17.78 ±0.12)ps- 

= (o.i54tEJ:EJ^,^) ps- 
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Using these inputs, together with the measured value of 
Ogj given in Eq. ([6|, we can extract the three theoretical 
parameters. We find a good fit, indicating that the mea- 
surements are compatible with one another, with result 



(8.889 ± 0.060)ps-i 
(0.112 ±0.040)ps-i 
-0.79 ±0.24. 
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From this one can see that the data prefers jMfjl to be 
close to the SM value, while both \T\2\ and 0'' differ 
from the values given in Eq. by about l.ba and 3ct 
respectively. This is in agreement with the result of [11] , 
which also found that a good fit to the data requires a 
non-zero phase as well as a value of \T\2\ higher than 
what is predicted in [S] . This is also compatible with the 
observation made in 1121, which found that new physics 
that only adds a relative phase 0* is unable to explain 
the central value of the semileptonic CP asymmetry. If 
we were to assume no new physics in the system, we 
would find the same value for iMfjl, but \T\2\ — (0.131 ± 
0.41)ps-i and 0" = -0.88 ± 0.24. 

Given this result, one might naturally be inclined to 
add new physics to [H].^ In the remainder of this 
paper we will study the constraints on NP contributions 
to from data on the decays of B mesons. The con- 
straints we derive are in general not sensitive to 0(1) 



factors neglected in our calculations. However, it is pos- 
sible that large numeric factors could relax or avoid some 
constraints. 

Any operator of the form 6si?, with R being any fla- 
vor neutral set of fields with total mass below tub^ can 
contribute to r^2- In order to conserve energy and mo- 
mentum, R needs to contain at least two fields. We first 
consider operators which only contain light fields present 
in the Standard Model, but comment on the possibility 
of introducing new light fields towards the end of the 
paper. The lowest dimensional operators possible have 
dimension six 



(13) 



where ip denotes any light Standard Model fermion. It 
is also possible to add a pair of operators, bsipiipj and 
bsipjipi, such that the combination is flavor neutral. A list 
of the possible operators is shown in Table |l] The physics 
of B decays is described by the electroweak Hamiltonian, 
which is conventionally written in the form 



H ^ 



(14) 



Characterizing the scale of new physics by A^p , we write 
the coefficients of the new operators as 



(15) 



Allowed operators 








Constr r 




Constr r 


bsuu 


K+H-, K+Tv° 


bduu 




hsdd 


K%+, K+7T° 


bddd 




bscc 




bdcc 




bsss 




bdss 




hsee 




bdee 


(tt, p)e+e^ 






bdfip. 


(7r,p)/i+^~ 


bSTT 




bdfr 


T+r- 


bsuv 




bdvi/ 


{■K,p)vV 


bssd 




bdsd 




bsds 


K°K°, K+K" 


bdds 




bscu 


DtTT-, K'^D'' 


bdcu 


D+Tv- 


bsuc 


D-K+, D°K+ 


bduc 





^ Note that the Standard Model predicts S^,p to be very close to 

zero, giving another hint at physics beyond the Standard Model 

in the Bs system. 
^ A recent measurement by the Belle collaboration ^ finds a 

value for AFs that is consistent writh the Standard Model values 

of irfjl and </-^ 

^ For previous attempts to explain the CP asymmetry by new 
physics contributions to Mj'j , see 1121 115H17| . 



TABLE I: Possible operators of the form bqip^p, with 1/) being 
an SM fermion. In the second column we show some decays 
that can be used to constrain each operator. The next two 
columns show the same for operators in the Bd system, which 
are required to keep the Bd lifetime in agreement with the Bs 
lifetime. 

The contribution of an operator O^p to rj2 can be 
evaluated by performing an OPE. Comparing the result 
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with the dominant contribution in the SM, arising from 
the operator bscc with Wilson coefficient C Vet, we 
find 



IpNPI 
|J- 12 I 
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where we have neglected phase space factors. A new bscc 
operator that can interfere with the SM operator is an ex- 
ception, which we discuss in more detail below. In order 
for the contribution of new physics r^2^ to compete with 
the Standard Model contributions, the Wilson cocflicicnt 
of this new operator needs to satisfy 



NP 
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where A is the Cabibbo angle A ^ 0.2. This is satis- 
g^pm-^^/X. Note that we have neglected 



fied if A 
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< 



numerical factors that arise from the contractions over 
Dirac and color indices. While these can be substantial 
for certain operators, our conclusions are in general not 
affected by these factors. 

It is important to note that the operator O^p ^^^^ con- 
tribute to decays in addition to Bg decays through 
the parton- level process b — >■ s'ipip. While this decay is 
phase space suppressed, it is enhanced by two powers of 
n^B^/ and is in general the dominant contribution to 
the Bs width from this new operator. We can estimate 
its effect on the iJ^ and Bs widths, including phase space 
factors, by writing 
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NP 
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X fd,sim^,/mb) . 
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Using the phase space factors given in [TH], an estimate 
of the function f{m^/mb) is shown in Fig. [ij The dif- 
ference between the two functions is dominated by the 
annihilation contribution of O^p ' which only affects the 
Bg system. The exception is a new bsdd operator, which 
contributes roughly equally to both Bd and Bg decay. 




0.1 0.2 0.3 0.4 0.5 mi, 

FIG. 1: The functions fd^s{m^/mt) 

One can easily see that for light fields the contribu- 
tion to the total lifetime of the Bd and Bg mesons can be 



as large as 50% or more. Given that these lifetimes have 
been measured to 1% and 2% accuracy, respectively, one 
might conclude that NP contributions to Ffj are com- 
pletely ruled out. However, our ability to predict these 
lifetimes accurately is plagued by large non-perturbative 
effects, and the resulting theoretical uncertainties could 
be as large as 10-20%. Therefore, we will attempt to 
find alternate bounds on NP contributions to these oper- 
ators, especially in the cases where ip = cott, for which 
the contributions are suppressed by phase space factors. 

To constrain these operators, we need to consider how 
they contribute to observable decays of the B mesons. 
Any non-leptonic operator of the form 651^293 will con- 
tribute to non-leptonic B decay B — )■ M1M2, where the 
fiavors of Mi and M2 depend on the flavors of the quarks 
Qi. While non-perturbative effects make it difficult to 
predict the precise rate for non-leptonic decays, factor- 
ization theorems exist at leading order in l/m>. |19l 120) . 
This allows us to estimate the decay rates as 



Bt{B M1M2) 



tbGI\C\ 



2 flimlFB^M 
327r 



(19) 



where C denotes the Wilson coefficient of the given oper- 
ator. Here we have used the scaling C ^ A^ and the rough 
estimates /m ~ 0.15 GeV and Fb^m ~ 0.3 to obtain a 
numerical value for the decay rate. Leptonic operators 
will contribute to decays of the form B M£^l^ , with 
a branching ratio estimated to be 21J 



Br(S ^ Mit+t 



'^sGf|C| ino_.^ PS{mi/mb) 



1927r3 
0.02PS{me/mb) , 



(20) 



One finds for the phase space factor PS(0) = 1 and 
PS(mT-/mh) = 0.05. They also contribute to the anni- 
hilation decay [2T] 



Br(B ^ e+r 



rBGl\Cf^4^Him,/mb) 

O/TT 

0.3 H{me/mb) , (21) 



where Hirni/nib) is a helicity suppression factor that is 
I ml if the the decay is helicity suppressed and unity 
otherwise, and we have used fs ~ 0.24 GeV for the nu- 
merical estimate. 

Finally, operators of the form bsijiil) can also con- 
tribute to the decay B Xg^ by mixing into the op- 
erator O7, which mediates this decay. If the NP oper- 
ator is an (axial)-vector current, this mixing occurs at 
0{as) {0{a)), but is enhanced by large logarithms of 
nT-b/rrit^w- Numerically, a non-leptonic (axial)-vector op- 
erator with a coefficient of order GpVcb gives a contribu- 
tion to S — >■ Xgj of the same order as the SM contribu- 
tion, while a leptonic operator would be suppressed by a 
factor of a/ag. If there is a (pseudo)scalar or tensor con- 
tribution, the mixing can occur at leading order and is 
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still enhanced by large logarithms [52] . This gives a con- 
tribution which is enhanced by 47r/as compared to the 
SM contribution"*. Of course, one can make the helicity 
structure of the NP operator such that mixing into O7 
is forbidden, but such operators will mix instead into an 
O7, which still contributes to i? — > Xsj- The resulting 
contribution to the branching ratio can be estimated as 



Br(B^M7)NP 
Br(B^M7)SM 



Vcb 



NP 
Vcb 



where C^p is the Wilson coefficient of the operator that 
mixes with the SM O7 and C^p corresponds to the op- 
erator that mixes with O'^. The variable r is 1 for 
non-leptonic (axial)-vector operators, a /as for leptonic 
(axial)-vector operators, and -iT:/ as for (pseudo)scalar 
and tensor operators. 

Given these results, one can immediately rule out the 
operators with ip = e, 11, since the observed branching 
ratios of B ^ K'-*h+l- are measured to be 0(10 
not O{l0^^) as the presence of new operators would pre- 
dict. The same is true for -0 = due to the limit on the 
branching ratio B — )■ K^*^Di>. Thus, of the leptonic oper- 
ators, only ip — T is allowed. The non-leptonic operators 
with light quarks are all ruled out by the absence of any 
2-body non-leptonic B decays to light mesons (such as K 
and tt) at the 10~'^ level. This also rules out the combi- 
nation of bscd and bsdc. The operator bscc |23j . however, 
cannot be excluded by this argument, since there is an 
SM contribution to this operator at the same level and 
the presence of non-perturbative effects makes a detailed 
comparison difficult. 

Both the remaining cases (^ = r, c) can be constrained 
by considering their contributions to the decay B Xsj- 
Given that the measured value f4j of Br™(_B Xsj) = 
3.52 ± 0.25 X 10~^ is consistent with the theoretical pre- 
diction [21] of Br"^(S ^ Xs-f) = 3.15 ± 0.23 x lO^*, 
only an 0(10%) correction can be accommodated. This 
eliminates any operator of the form {bs)(ijj'ip)s.p,T, due 
to the factor Air /as in Eq. (22). Note that this includes 
the operator discussed in (14j . 

The operators (bs){fT)v,A can not be constrained 
because, as discussed, the mixing only occurs at 
one loop and is suppressed by a /as- The operator 
(bs)v~A{cc)v±Ai which mixes with the operator O7, 
can be eliminated, since its contribution to the decay 



B Xsj is of order \r^{\ / 



l-nSMI 
P- 12 



An operator that 



with helicity structure {bs)v+A{cc)v-A has this property 
and can therefore contribute significantly to the lifetime 
difference in the Bs system. 

We have shown that there are only two possible SM 
operators that can give rise to an 0(1) change in Ffj- 
The first is (bs){fT)v,A- This operator can be constrained 
by both B — > i^'^*V+T~ and Bs — ?> t+t~; however, due 
to the difficulty in detecting r's, there is currently no 
bound on either decay. We therefore find that bsfr can 
contribute significantly to to Ffj- The second possible 
operator is of the form (bs)v+Aicc)v-A- 

Note that both of these operators would give rise to an 
order 10% contribution to the total lifetime of the Bg me- 
son, if we require that they contribute an 0(1) amount to 
rf2- As discussed above, however, this does not contra- 
dict the precise measurement of the Bs lifetime, due to 
the large theoretical uncertainties when trying to predict 
this quantity. On the other hand, the ratio tb, / tb^ is un- 
der much better theoretical control. This is because the 
unknown nonperturbative effects largely cancel in the ra- 
tio, such that it can be predicted with high accuracy [IS] 



r{Bs) 



= 1 ± 0(1%) . 



(23) 



An operator that gives an 0(1) contribution to would 
give rise to large lifetime difference 1 — t{Bs) /T{Bd) — 
0(10%), much larger than the theoretical uncertainty in 
this quantity. Unfortunately, the experimental uncer- 
tainties in the ratio [10] 



r{Bs) 



0.965 ±0.017 



(24) 



are somewhat larger than our theoretical knowledge. 
While one can rule out a 10% effect, a 5% contribution 
is still allowed. In fact, the current measurement seems 
to indicate a 2a deviation in the lifetime ratio. If a sig- 
nificant difference from unity of this lifetime ratio could 
be established, it would be another hint at new physics 
contributing to Ffj- A more precise measurement of this 
quantity is therefore of great importance. 

What would one conclude if a new measurement of this 
lifetime ratio does not allow for a large deviation from 
unity? Since the operators discussed above reduce the 
ratio Bs lifetime relative to the Bd lifetime, one would 
be forced to add new operators of the form 



Qd 



NP 



mixes only with O7, on the other hand, contributes only 
quadratically to S -> Xs^f, and therefore C^p/Kf, 0.3 ^j^i^ Wilson coefficient 
would still be allowed. For this to lead to a sizable effect 
in F^2j the operator has to interfere with the SM operator 
{bc)v-A{cc)y-A in its contribution to Ff2- An operator 



6F(iV'iF'02 



NP 



NP J 



(25) 



(26) 



We thank Uli Haisch for discussions on this point. 



to make up for this difference. In general, we do not 
need the fields in O^p to be the same as in O^P' which 
allows ip ipi ip2- Note that as long as ipi 7^ ■02, 
such operators would not contribute to Tf2- While one 
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could potentially add several new operators to the B^, 
sector with smaller Wilson coefficients to compensate for 
Onp j we will assume that only one new operator is added. 
Our arguments can be easily generalized to the case of 
multiple operators. Adding operators to the Bd system 
would also change the ratio r {B^) /r {Bd), which is also 
well understood theoretically [57], [3S]. However, we will 
not explore new operators in the i?„ system here. 

The possible operators are shown on the right hand 
side of Table |l] As before, the operators bdll for I — e,^ 
are ruled out by limits on the decay B — 7r/"*"Z^, and 
the operator with ijj — hy the limit on B ^ ttDi'. The 
operator bdfr is in this case excluded by the experimen- 
tal limit on the decay B — T+r~. Operators involving 
only light quarks are also ruled out, again due to the ab- 
sence of 2-body non-leptonic B decays to light mesons at 
the 10"'^ level. The operator bdcu (which has the wrong 
sign charm) is also excluded by this argument, while the 
operator 6dcc is ruled out by its large contribution to 
B^(p,a.")7. 

This leaves only one operator one could add to the 
Bd system, namely bduc. A new physics operator with 
C^p ^ would contribute at the same order as the 
SM operator. While this would increase the rate of non- 
leptonic B decays of the form B — Dir, the presence of 
non-perturbative effects makes it difficult to rule out this 
operator conclusively. Note that this operator makes it 
possible to keep the lifetime ratio between Bs and Bd 
mesons the same, as long as its coefficient is tuned suffi- 
ciently Given that the operator bduc appears to be com- 
pletely unrelated to either the bsfr or (bs)v+Aicc)v-A 
operators allowed in the Bg system, such a NP scenario 
seems very contrived. 

Having discussed in detail the effect of dimension six 
operators containing SM fields, we want to briefly com- 
ment on other possibilities for NP, while stressing that 
we can not rigorously discuss all possible extensions. The 
simplest extension would be to allow for new light degrees 
of freedom (denoted by (f> and for spin and spin 1/2), 
that are SM singlets (or electrically neutral, but charged 
under SU(2)l). Operators of this form include b{s,d)(j)<j) 
and b{s,d)'iptp, b{s,d)(j)F and b[s,d)(pG where F and G 
denote the photon and gluon field strength, respectively. 
The first two operators are ruled out by the experimen- 
tal bounds on B ^ {K'^*\uj, p)Dv. The operator bd4>F is 
ruled out by the absence of the decay B — vv^, but the 
other operators are not ruled out, at least not without a 
more careful analysis. 

As another possibility, one might consider higher di- 
mensional operators containing SM fields. At dimension 
seven, a derivative could be added to any of the oper- 
ators we have already discussed. All of the previously 
mentioned constraints would still need to be considered; 
in addition, it is unlikely that dimension seven opera- 
tors could contribute at the same level as the dimen- 
sion six operators. However, three new operators ap- 



pear: b{s,d)FF, b{s,d)GG and b{s,d)FG, where F and 
G are the photon and gluon field strengths respectively. 
The operators with two photons are ruled out by the de- 
cays Bd^s 77 and B Xd^sl- The operators with 
two gluons would generate four quark operators of the 
form b{s,d)qq of roughly the same size as the previously 
studied dimension six operators and are thus excluded. 
Finally, the operators with one photon and one gluon 
also contribute to B Xd.sJ, and are therefore ruled 
out. One could continue to add higher dimensional SM 
operators, but given that we need a fairly large contribu- 
tion to rf2, it is unlikely that this could be done without 
lowering Anp into a region that is excluded by collider 
experiments. 

As a last possibility, we would like to comment on op- 
erators that contain one or more new light degrees of 
freedom, but which eventually will decay to SM parti- 
cles. Assuming that light degrees of freedom that are 
charged under SU(3)c or have 0{e) electric charge are 
already ruled out by collider data, we are left with fields 
that couple to 6(s, d) through higher dimensional oper- 
ators. They also have to couple to the Standard Model 
in such a way that they predominantly decay to particles 
charged under the SM before exiting the detector, with- 
out having any 0(1) SM charges. Finally, these operators 
are constrained to contribute to Ts^d at no more than the 
20% level and maintain the measured ratio tbJtb^- It 
might be an interesting exercise to see if such a scenario 
can be realized by clever model building. 

In conclusion, we have analyzed if current data al- 
lows significant new physics contributions to Ffj, as re- 
cent measurements might suggest. We have shown that 
there are two dimension six operators in the Standard 
Model which are still allowed given current experimental 
constraints. More detailed measurements of the decays 
B — > XgT^T^ and Bs — > t+t^ would establish the pres- 
ence of one of these operators or rule it out. Both of these 
operators would contribute differently to the B^ and the 
Bd lifetime. This ratio is predicted very accurately in the 
SM, and an observed value in agreement with the Stan- 
dard Model prediction (and similarly small uncertainties) 
would require additional physics in the Bd system in a 
way that seems unrelated to the new physics allowed in 
the Bg system. Therefore, a more precise measurement 
of T{Bs)/T{Bd), which is possible with currently avail- 
able data, could either support or strongly constrain the 
existence of new physics in Ff 2 . 
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